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Introduction: On Earth, mass wasting events such 
as rock falls and landslides are well known conse- 
quences of seismic activity. Through a variety of re- 
mote sensing techniques, tectonic faults have been pos- 
itively identified on all four of the inner planets, 
Earth’s Moon, several outer planet satellites, and aster- 
oids [1]. High-resolution imaging has furthermore ena- 
bled positive identification of mass wasting events on 
many of these bodies. On Mars, it has been suggested 
that fallen boulders may be indicative of pale- 
omarsquakes [2]. On the Moon, meteor impacts and 
moonquakes have likewise been suggested as potential 
triggering mechanisms for mass wasting [3]. Indeed, 
we know from the Apollo era that the Moon experienc- 
es a wide variety of seismicity [4] . 

Seismicity estimates play an important role in creat- 
ing regional geological characterizations, which are 
useful not only for understanding a planet’s formation 
and evolution, but also of key importance to site selec- 
tion for landed missions. Here we investigate the re- 
gional effects of seismicity in planetary environments 
with the goal of determining whether surface features 
such as landslides and boulder trails on the Moon, 
Mars, and Mercury could be triggered by fault motion 
(Fig. 1). We attempt to quantify the amount of near- 
source ground shaking necessary to mobilize the mate- 
rial observed in various instances of mass wasting. 

Lobate scarps: Lobate scarps, the typical surface 
expressions of thrust faults resulting from tectonic 
compression, are widely observed on the Moon, Mars, 
and Mercury (Fig. 2). Compared to other types of tec- 
tonic faults, surface-cutting thrust faults require the 
largest amount of stress to form and/or slip, and thus 
they could be expected to result in large quakes. While 
normal faults, graben, and wrinkle ridges may be more 
abundant on Mars, the Moon, and Mercury respective- 
ly, these structures would create smaller theoretical 
maximum quakes than lobate scarp thrust faults. Thus, 
we optimize our chances of finding mass wasting asso- 
ciated with faults by studying lobate scarps. 

Methodology: Given an observed fault, we first 
focus on calculating the theoretical maximum quake 
that could occur as a result of slip there, and then de- 
termine the resulting effects on the surrounding surface 
morphology. The expected damage area indicated by 
seismic wavefield modeling can be compared to 
mapped imagery to determine the likelihood of a quake 
having triggered a mass wasting event. 



Fig. 1: (left) Landslide deposits (granular flow) on an inte- 
rior slope of Marius crater on the Moon (11.9°N, -50.8°E). 
(right) Boulder tracks emanating from a crater rim alcove 
on Mars (-9.515°N, 16.433°E). A 74-km compressional 
fault in the Arabia-Sabaea Terra is located < 100km away. 



Fig. 2: Examples of lobate scarps on the Moon (left), Mars 
(center), and Mercury (right). Moon: Ever shed SI (center 
lat/lon 33°N/197.1°E), Mars: Utopia Planitia#s 1801, 1802, 
1804 (center lat/lon 52.9°N/119.2°E), Mercury: Beagle 
Rupes (center lat/lon -3.5°N/100.7°E). 


Theoretical maximum quake. Following the meth- 
od outlined in [5], the theoretical maximum quake 
magnitude is derived from basic fault properties. These 
are either estimated from imagery or derived from la- 
boratory rock experiments or elastic dislocation mod- 
els, and include the length (L), dip angle (8), depth of 
faulting (T), and fault width (w) (Fig. 3). Fault dis- 
placement (D) is calculated using displacement-length 
scaling such that D = yL, where y is determined by 
rock type and tectonic setting [6]. We note that subsur- 
face fault geometry and mechanical properties of 
planetary lithospheres and regoliths are incompletely 
understood, and thus represent potential sources of 
error in the maximum quake calculation. To incorpo- 
rate this error, we investigate ranges in fault parame- 
ters, placing upper and lower bounds on our maximum 



Fig. 3: Schematic 3-D (left) and cross-sectional (right) 
views showing the fault parameters: displacement (D), dip 
angle (8), vertical relief, depth of faulting (T), and fault 
width (w) for the thrust fault underlying a lobate scarp. 


quake calculations rather than estimating discrete val- 
ues. 

The best measure of the size of a planetquake is its 
seismic moment, Mo. It is calculated by multiplying the 
shear modulus of the ruptured rock (G) by the area of 
the ruptured portion of the fault (A) and the average 
displacement (D) produced during the quake, such that 
Mo = GAD = G(Lw)(yL). The seismic moment repre- 
sents the total energy consumed in producing dis- 
placement on a fault, regardless of the local strain rate 
or fault formation mechanism. 

Seismic wavefield modeling. In order to determine 
the dimensions of an area affected by seismic shaking, 
we model the ground motion resulting from the theoret- 
ical maximum quake along a given fault (Fig. 4). Fol- 
lowing the method of [7], we use the the Serpentine 
Wave Propagation Program (WPP), a numerical code 
for simulating seismic wave propagation through arbi- 
trary elastic and anelastic media in a 3-D model space. 
The initial model of a given fault includes regional 3-D 
topography derived from digital elevation models, and 
the planet’s relevant background 1-D velocity model. 

We note that the modeled peak ground motion is 
not as strongly dependent upon the choice of back- 
ground velocity model, but more so upon the scattering 
and attenuation properties of the shallowest materials 
in the model. Synthetic seismograms for the Moon 
most reasonably approximating those recorded by the 
Apollo seismometers are acquired for a 1 km thick, 
highly scattering layer as the topmost layer in the mod- 
el. Similar highly fractionated layers are expected on 
Mars and Mercury, and we approximate their velocities 
using the physical properties of a basaltic crust for each 
of body. 

Testing: Peak vertical ground velocity (as a proxy 
for displacement) typically occurs within a few kilome- 
ters of the main shoch and drops off rapidly from there. 
This implies that we should expect most of the land- 
slides and other mass wasting phenomena to occur in 
the immediate vicinity of the fault. However, this result 
may depend on regional effects such as surface slope 
and megaregolith thickness. For example, a thicker 
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Figure 6: Predicted maximum vertical ground motion in the 
vicinity of the Evershed lobate scarp on the Moon using the 
3-D WPP code. Left: Surface topography input into the 
simulation, from the Lunar Orbiter Laser Altimeter experi- 
ment. The Evershed scarp is centered in the image 
(162.9°W, 33°N, see Lig. 2). Right: Ground motion for a 
Mw=7.8 quake on a subjacent reverse fault, with a 2.25 km 
depth of faulting. The surface trace of the scarp is indicated 
by the red line. A random distribution of heterogeneity of 
25% in shear and compressional wave velocity with 100 km 
scale length scatterers is placed in the lunar megaregolith to 
simulate the seismic scattering typically present in lunar 
seismograms. Peak ground velocity is measured for the first 
1000 seconds of the seismic trace. 


megaregolith (as might be expected in the vicinity of 
craters) would tend to focus shaking in some of the 
crater basins. The presence of sediments also enhances 
seismic shaking; this could be a relevant scenario for 
Martian craters that may have been lakes some time in 
the past. 

We will compare the observed extent of mass wast- 
ing in the vicinity of a fault to the modeled event mag- 
nitude and peak ground motion in order to establish a 
method to translate quake parameters into mass wast- 
ing estimates. This has been perfomed for terrestrial 
examples focused on determining landslide area and 
density over time in seismically active regions, as well 
as using the presence or absence of precariously 
perched boulders as indicators of the regional seismic 
shaking. The latter example has also been performed 
on Mars, where both boulder size and boulder trail 
density were found to peak close to the center of a fault 
system and decrease linearly along strike [2]. We ex- 
pect to find systematic variations in fit parameter esti- 
mates for each body, reflecting different gravitational 
strengths, regolith cohesion properties, and other geo- 
logic settings local to each body/study region. 
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